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ABSTRACT 

We investigate the dynamical response, in terms of disc size and rotation velocity, to mass 
loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin 
density profile embedded in a spherical dark matter halo having a density profile proposed 
by Navarro, Frenk & White is considered. For a purpose of comparison, we also consider 
the homogeneous and r _1 profiles for dark matter in a truncated spherical halo. Assuming 
for simplicity that the dark matter distribution is not affected by mass loss from discs and 
the change of baryonic disc matter distribution is homologous, we evaluate the effects of 
dynamical response in the resulting discs. We found that the dynamical response only for an 
adiabatic approximation of mass loss can simultaneously account for the rotation velocity and 
disc size as observed particularly in dwarf spiral galaxies, thus reproducing the Tully-Fisher 
relation and the size versus magnitude relation over the full range of magnitude. Furthermore, 
we found that the mean specific angular momentum in discs after the mass loss becomes larger 
than that before the mass loss, suggesting that the mass loss would occur most likely from the 
central disc region where the specific angular momentum is low. 

Key words: galaxies: spiral, disc - galaxies: dwarf - galaxies: evolution - galaxies: formation 
- galaxies: haloes - large-scale structure of the universe 



1 INTRODUCTION 

It is well-known that the luminosity L of luminous spiral galax- 
ies correlates tightly with their rotatio n velocity V TO t, as called 
the Tully-Fisher relation (hereafter TFR. lTullv & Fished 1977b . The 
TFR is usually described as a power-law scaling relation L oc V^ ot , 
and the rms scatter of the TFR is smaller for optical passbands 
of longer wavelengths. In particular, the I-band TFR is signifi- 
cantly tight and its scatter is approximately as small as 0.4 mag 
(Willick et al. 1997; Kannappan, Fabricant & Franx 2002), with a 
powe r-law index 7 ~ 3 — 3.5 (e.g. lPierce &Tully 1988;Sakaietal. 
2000). Light at longer wavelengths from galaxies traces old late- 
type stars, unaffected by sporadic formation of young stars, and 
represents the mass of luminous matter. Thus, the tightness of the 
I-band TFR suggests the existence of some driving mechanism that 
depends on mass of galaxies in processes of their formation. 

The current standard model of cosmology is the cold dark mat- 
ter (CDM) model. Since the initial spectrum of density fluctuations 
has larger amplitudes for smaller scales in this model, the scenario 
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of structure formation is hierarchical clustering, in the sense that 
smaller dark haloes cluster to form larger dark haloes hierarchi- 
cally, thereby creating a scaling rel ation between mass M and cir- 
cular velocity V^i rc of dark haloes (Faber 1982; Blumenthal et al. 
1 1984) . High-resolution iV-body simulations, based on a power 
spectrum P(k) oc k n with an index n = —2.5 appropriate for 
galaxy scales or reciprocal of wavenumber have provided 
M oc V^ ir , reminiscent of the TFR if a constant mass-t o-light 
M/L ratio is assumed (e.g. iNavarro. Frenk & Whitd [i997). Con- 
sequently, the CDM model naturally involves a TFR-like scaling 
relation as observed for luminous spiral galaxies, because a con- 
stant M/L i s indeed known to be a g ood approximation to such 
galaxies (e.g. lFaber & Gallagheal 19791) . 

When studying the formation of galaxies originated from the 
growth of density fluctuations in the early universe, however, it is 
also necessary to take into account the effects of star formation 
and supernova (SN) explosions in individual galaxies. In particu- 
lar, dwarf galaxies have shallow gravitational potential wells, so 
that energy feedback from SN explosions significantly affects their 
evolution and hence scaling relations that follow. In a framework 
of monolithic collapse scenario of galaxy formation, several au- 



2 H. Koyama, M. Nagashima, T. Kakehata, and Y. Yoshii 



thors studied scaling relations for dwarf galaxies in two extreme 
cases, i.e., the dark matter dominat ed case in which d ark halo 
dominates the gravitational potential dDekel & Silklll986h . and the 
baryon dominated case of self -gravitating galaxies without dark 
matter dYoshii & Arimotolll987l) . 

ISteinmetz & Navarro! (l999) derived, for the first time by nu- 
merical simulations, the TFR for luminous spiral galaxies in the 
CDM model, using a method of TV-body and smoothed particle hy- 
drodynamics (SPH) combined with star formation and SN feedback 
explicitly. They also showed that the star formation rate is regulated 
for M/L to be almost constant. However, the zero-point of the TFR 
turned out to differ from that observed, so that the simulated spi- 
ral galaxies were too faint to be consistent with observations at 
the same circular velocity. It was therefore pointed out that pro- 
cesses of star formation and SN feedback i n current iV-body/SPH 
simulations have to be improved (see also iGovernato et alj|2007l ; 
IPortinari & Sommer-Larsenll2007l) . 

Semi-analytic (SA) models of galaxy formation have been 
invented as a complementary approach to numerical simula- 
tions, in which complex processes such as star formation and 
SN feedback are simply modeled on galaxy scales. SA mod- 
els have indeed succeeded in reproducing many observational re- 
sults, for example, the luminosity function of galaxies as well 
as the TFR, by sup pressing the formation of dwarf galaxies ow- 
ing to SN feedback dSomerville & Primackll 19991: ICole et all200(j| ; 
Nagas hima & Yoshiil |2004 iDe Lucia. Kauffman n & White! |2004 



Kang et all2005l ; ICroton et al.ll2006h . 



Recent progress in observational techniques has enabled to see 
much more properties of distant galaxies and faint dwarf galaxies 
beyond the limit of magnitude attained so far. In contrast to the suc- 
cess of SA models in accounting for luminous galaxies, they came 
to face a serious discrepancy between predicted and observed dy- 
namical properties of dwarf galaxi es, while reproducing their pho- 
tomet ric properties quite well fe.g. lCole et a l. 2000; van den Boschl 
|2002|) . 

Here we would like to stress that the dynamical response to 
gas removal induced by SN explosions, which has been overlooked 
in most SA models, is unavoida ble in dwarf, less massive galax- 
ies. iNagashima & Yoshiil J2003h formulated how the structure of 
spherical galaxies responds to SN -induced gas removal in a gravi- 
tational potential of dark halo, and Nagashima & YostriJ d2004l) in- 
corporated this effect into a SA model to show that many prop- 
erties of elliptical gala xies, including the Faber-Jackson relation 
dFaber & Jackson! 1 19761) . can be reproduced to faint magnitudes 
of dwarf ellipticals. A more sophisticated SA model associated 
with high-resolution TV-body cosmological simulations also pro- 
vides a good ag reement between pred icted and observed proper- 
ties of galaxies ([Nagashima et al. 2005). Accordingly, as a natural 
extension of Nagashima & Yoshii's trial, it is worth investigating 
whether such a dynamical effect could also reproduce the observed 
TFR down to faint magnitudes in spiral galaxies. 



In this paper, using the Kuzmin disc jKuzmin|[T952lll956h 
a galaxy disc embedded in a spherical dark halo, we present in Sec- 
tion 2 the general formulation for the dynamical response in size 
and rotation velocity of discs. In Section 3 we examine such effect 
for several choices of density distribution of dark halo. In Section 4 
we derive the resulting TFRs and disc size versus magnitude rela- 
tions. In Section 5 we summarize the results of this paper. Detailed 
derivations of analytic expressions are given in Appendix. 



2 DYNAMICAL RESPONSE OF GALACTIC DISK IN 
DARK MATTER HALO 

In this section, we formulate the dynamical response accompanied 
by SN feedback to a spiral galaxy consisting of baryons and dark 
matter. We consider a thin baryonic disc or the Kuzmin disk embed- 
ded in an non-rotating spherical dark halo whose density profile is 
either NFW, homogeneous or r~ . We assume that the thin disc 
is axisymmetrical and rotation-dominated with negligible velocity 
dispersions. Note that the density profile of the Kuzmin disc is more 
or less similar to the exponential disc except for the central region. 

The gas of low angular momentum in an extended halo col- 
lapses towards the galaxy centre on a dynamical timescale and dis- 
sipates the energy to cause a burst of star formation in the central 
region. On the other hand, the gas of higher angular momentum 
gradually falls onto a disc plane on a longer timescale and settle on 
a circular annulus of the disc at the radius according to the value 
of angular momentum of the gas. This fallen gas is converted into 
stars to form a stellar disc. 

If the kinetic energy, which is released by SNe following the 
star formation, exceeds the binding energy of the disc, the remain- 
ing gas is removed out of the disc. This disc dynamically recovers a 
final equilibrium. In general, the final state could be either a puffed- 
up disc where velocity dispersions dominate over the rotation, or a 
thin disc where the rotation remains to domina te with negligible 
velocity dispersions dBiermann & Shapiroll979T) . In this paper, it is 
enough to consider only the latter, because we are primarily inter- 
ested in scaling relations in spiral galaxies. 

First, we consider that the gas removal occurs after all the ma- 
terial falls in to form the disc. In this case, the disc achieves a cen- 
trifugal equilibrium with the mass Mi and the angular momentum 
Ji as an initial state. Since the gas removal accompanies the si- 
multaneous losses of mass and angular momentum, the disc should 
newly recover a centrifugal equilibrium with the mass M f and the 
angular momentum Jf as a final state, depending not only on the 
amount of such losses AM(Mi — Mf) and AJ(Ji — Jf), but also 
on whether the time scale of such losses is lon ger or shorter than 
the dynamical scale dBiermann & Shapiroll 19791) . 

On the other hand, the gas removal may well occur while the 
material is still falling into the disc. In this case, only a part of the 
material falls in to form the disc and almost no gas removal occurs 
afterwards. The disc would then be settled in centrifugal equilib- 
rium with the mass Mf and the angular momentum Jf throughout 
from the beginning. This situation is formally equivalent to the case 
such that the gas removal occurs on much shorter timescale com- 
pared to the dynamical timescale. Therefore, the related formula 
below in this section could also apply to the situation considered 
here. 



2.1 Formulation 

The kinetic energy T due to rotation, the gravitational self plus 
interaction potential energy W, and the total angular momentum J 
of the Kuzmin disc of baryons in a dark halo are given by 



T= ^M b V b 2 , 

w = - 1 7 gMl 

4 r b 

and 



J = gM b r b V b , 



„M b M d 
- G f(z), 



(1) 



(2) 



(3) 
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respectively, where G is the gravitational constant and g is the con- 
stant dependent on density distribution, M is the mass, V is the 
characteristic rotation velocity, and r is the characteristic radius. 
The subscripts "6" and "d" refer to the baryonic and dark matter 
components, respectively. The function f(z) has a form dependent 
on density distribution, and the argument z is a radius rt, normal- 
ized by r d or z = r b /r d . 

Hereafter, the quantities of baryons with the subscript "i" in 
place of "6" are designated as an initial state before the gas re- 
moval, and those with the subscript "/" as a final state after the gas 
removal. Note that the quantities of dark matter with the subscript 
"d" are assumed not to change throughout the gas removal. 

First, we describe the initial state expressing the total energy 
of the disc as 



1 1 /Vf 2 
Ei=Ti + Wi= ~MiV? - -G — - 

2 4 n 



G^Mlf(z, 
rd 



(4) 



provided that the virial equilibrium 2Ti = —Wi holds, and the 
angular momentum as 



Ji = gMinVi. 



(5) 



Next, assuming that the disc profile does not change just after 
the gas removal, we express the total energy as 



Ef = ^MfVi 



, 1 M[ 
4 n 



(6) 



Then, conserving this energy, the disc dynamically responds to- 
wards recovering a virial equilibrium as the final state: 



and 



\m s v] 



1 Mt 
4 r f 



^MfMd ., . 
G ? — I /(*/)> 



J f = gM f r f V s . 



(7) 



(8) 



We here consider an extreme approximation of impulsive mass 
loss such that the gas removal is faster than the dynamical relax- 
ation. From equations {4j> and Q, we express the final quantities in 
terms of the initial quantities as 



mj_ 

rm 

and 

Yl 

Vi 



l + 4(z l /m l )[f(z f )-f( Zi )} 
2 - Zi/zf 



m-f/zf + 4f(z f ) 



i/zi + 4/(zi) 



1/2 



(9) 



(10) 



where Zij = Tij/rd andmij = Mij/M d . 

We consider another extreme approximation of adiabatic mass 
loss such that the gas removal is slower than the dynamical relax- 
ation. The total change from the initial to final states is a sum of 
consecutive infinitesimal changes. Substituting m/ = rrn + dm 
and Zf = Zi + dz in equation {9]l and linearizing it, we obtain 



dm m , df(z) 

1T = + 4z ^ 

dz z dz 



(11) 



where the subscript "i" is omitted for simplicity. The first term 
comes from the self-gravitating baryons alone, and the second term 
from the gravitational interaction of baryons with dark halo. Solv- 
ing this differential equation, we obtain 

C 



+ q{z), 



(12) 



where 



dt 



(13) 



and C is an integration constant or an adiabatic invariant here. Then 
we express the final quantities in terms of the initial quantities for 
the adiabatic mass loss: 



mj_ 

mi 



C/z f + q(z f ) 



Zi_ q(z f ) - Ziq{zi)/zj 
Zf mi 



(14) 



together with the same form of equation 4 1 Ob for the rotation veloc- 
ity. 

It is apparent from equations ((9), d 1 Ot . and dl4t that the ra- 
tios of dynamical quantities r//n(= Z/) and Vf /Vi are obtained 
when the ratios of m,(= Mi/Md), 2<(= ri/ra), and m//mi(= 
Mf I Mi) are given. In other words, given m, and Zi, the ratios of 
r//r; and Vf /Vi are formally written as a function of M / jMc 



^ = H{mi, Zi ;Mf/Mi), 
n 

and 

Vf 

-f- = I(m l ,z l ;Mf/Mi). 



(15) 



(16) 



For the impulsive and adiabatic approximations of mass loss, 
the change of specific angular momentum is given by 



(J/M) f _ TfVf 



(17) 



{J/M)i rM ' 

from equations l[5} and l[8}. When the gas removal occurs from the 
central region of the disc where the mass is concentrated but the 
rotation velocity is small, (J/M) / /(J/M)i should exceed unity. 
On the other hand, when the gas removal occurs from the outer 
region of the disc where the mass is small but the rotation velocity 
is large, (J/M) / /(J/M)i is less than unity. Therefore we use the 
value of (J/M) f / (J/M)i to discuss where in the disc the site of 
efficient gas removal is located. 

2.2 Limiting cases 

In this subsection, before studying more general gas removal, we 
examine two limiting cases of Mi/Md S> 1 and Mi/Md <SC 1. 

In the baryon dominated case of Mi/Md 3> 1, we have for 
the impulsive mass loss 



11 

r, 

Yl 

v 



(J/M) f 
{J/M), 



MfJMj 
2M f /Mi - 



M f 

> L — i 

Mi L > 

Mf/Mj 
2M f /Mi - 



,Mf 

'Mi 



-1 > 1, 



(18) 



and for the adiabatic mass loss 



12 

n 

Yl 

v 

(J/M) f 
{J/M)i 



M_ 
Mf' 
Mf 
~~M~i'' 



= 1. 



(19) 



The fact that the value of (J/M) / /(J/M)i is always larger than 
unity for the impulsive mass loss suggests that the gas removal 
should preferentially occurs from the central region of the disc, 
whereas the fact of (J/M) / /(J/M)i = 1 for the adiabatic mass 
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loss yields that the same fraction of mass is lost across the disc over 
its entire radial range. 

In the dark matter dominated case of Mb/Md <C 1, the disc 
size and rotation velocity of the initial disc hardly change during 
the gas removal, because the dominant gravitational potential of 
dark halo allows almost no response in structure and dynamics of 
the disc. 



3 DYNAMICAL RESPONSE ON DISK SIZE AND 
ROTATION VELOCITY 

In this section, based on the analytic expression of f(z) for the 
Kuzmin disc embedded in a dark halo with the NFW, homoge- 
neous, and 1/r density distributions, we calculate the strength of 
dynamical response for rf/n, Vf/Vi, and (J/M) / /(J/M)i as a 
function of Mf /Mi for several combinations of rrii = 0.05 — 0.2 
and Zi = 0.05 - 0.2. 

The range of mi explored is based on the following consider- 
ation. We assume that all the baryonic matter is distributed in the 
same way as the dark matter initially and falls in to make up the 
disc. Then the upper bound of rrn is set by the current estimate of 
density ratio pb/pd — 0.2 in the universe. Considering a possibil- 
ity that not all of the baryonic matter falls into the initial disc, we 
rather arbitrarily adopt a range of rrii = 0.05 — 0.2 for the initial 
mass ratio. 

The range of Zi = n/ Td is based on the following consider- 
ation. From the observed luminosity profile and rotation curve of 
dwarf spirals, the scale radius of dark halo of 10 9— 10 M@ is esti- 
mated as (1.5 — 4)/i _1 kpc, and the correspondin g rotation veloc- 
ity is estimated as 40 - 70kms _1 teurked[l995ri . where h is the 
Hubble constant Ho defined as h = i?o/100kms _1 Mpc _1 . In ad- 
dition, the radial exponential scalelen gth of such dwarf spirals is 
estimated as (0.5 - 2 )fe -1 kpc jPildis. Schombert & Edgerll 19971 : 
ISchombert et alj[l997h . thus giving the scale ratio of 0.3 — 0.5 as 
observed in the local universe. Regarding this ratio as obtained after 
the initial disc expands by a factor of several due to the gas removal 
(see below), we adopt a range of Zi = 0.05 — 0.2 for the initial scale 
ratio. 



3.1 Kuzmin disc in the NFW halo 

We examine the dynamical response to the Kuzmin disc embed- 
ded in the NFW halo. The surface mass density distribution of the 
Kuzmin disc is given by 



E(r) 



i r /<5$ 

AnG \ V Sz J z — o+o V Sz S z-,0-0 
M b r b 



2 7r(r 2 + 7-2)3/2 ' 
with the gravitational potential 

. GM b 
$(r, z) 



y/r 2 + (r b + \z\) 2 
As for the NFW halo, the density distribution is given by 



p(r) = p d c s 



with the gravitational potential 

*(r) = -4tvG^ In f 1 + -) , 
r \ r d J 



(20) 



(21) 



(22) 



(23) 



where c is a concentration parameter. Since the essential charac- 
teristic length of the NFW halo is r^/c instead of ra, we hereafter 
define z = r b /(r d /c). 

Figure Q] shows the dependence of dynamical response on rrii 
and Zi with c = 10 for the impulsive mass loss in the left column, 
and for the adiabatic mass loss in the right column. In each of the 
columns, the strength of dynamical response for r//rj, Vf /Vi, and 
(J/M) f / (J/M)i is shown in order from top to bottom panels as a 
function of Mf /Mi for several combinations of rrii and Zi. 

It is clear from this figure that the dynamical response 
monotonically increases as Mf /Mi decreases and asymptotically 
reaches a constant level at Mf/Mi ^ 0.1. Furthermore, given 
Mf /Mi, the dynamical response is stronger for increasing the ra- 
tio of surface mass densities of initial disc and dark halo projected 
onto the disc plane or rrii/z 2 . Such a tendency holds irrespective 
of either the impulsive or adiabatic mass loss. However, comparing 
the left and right columns, we find that the dynamical response is 
greatly suppressed for the adiabatic mass loss, in particular, even 
by a factor of ten or so within the range of rrii/z 2 and Mi/Ma 
considered. 

Figure [2] shows the dependence of dynamical response on c 
with rrii = 0.2 and Zi = 0.05. There is a tendency such that the 
dynamical response is stronger for increasing c, or equivalently de- 
creasing Ti/r d because Zi = r»/(rd/c) is fixed here. 

The specific angular momentum after the gas removal 
(J/M) f is always larger than that before the gas removal (J/M)i, 
which suggests that the gas removal occurs from the central region 
of the disc where specific angular momentum is small. This is the 
condition that the disc should remain thin and rotation dominated 
before and after the gas removal. 



3.2 Dependence on dark matter distribution 

We examine how the Kuzmin disc responds depending on the 
density distribution of dark halo. Here we consider the homoge- 
neous and 1/r density distributions in a spherical halo truncated at 

r = r d : 



P(r) = Pd0(r d - r), 



and 



p(r) = Pd — 0(r d - r), 
r 



(24) 



(25) 



respectively, where 6(x) is the Heaviside step function. 

The results for these two distributions are shown in Figures [3] 
and|4] respectively. We see from these figures that, given m, and Zi, 
the dynamical response for the homogeneous density distribution is 
stronger than that for the 1/r density distribution with which dark 
matter is more confined inside the characteristic radius of the initial 
disc Vi. 

We also see that, given rrii and Zi, the dynamical response for 
the NFW halo (Figure QJ is even stronger than those considered 
here. This is because, adopting c = 10 in the definition of Zi = 
n/(rd/c), the NFW halo is more extended beyond the truncation 
radius ra, thereby dark matter is less confined inside the initial disc 
radius n as far as the value of rrii is fixed. 

Furthermore, the dynamical response for the impulsive mass 
loss is much stronger than that for the adiabatic mass loss, which is 
also the case of the NFW halo (FiguresQ]and|2]l. 
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Figure 1. Dynamical response of the Kuzmin disc embedded in the NFW halo. The results shown are the dependence of dynamical response on rrii = Mi/Mj 
and Zi = Ti/(r^/c) with c = 10 for the impulsive mass loss in the left column, and for the adiabatic mass loss in the right column. In each of the columns, 
the strength of dynamical response for rj/ri (top panel), Vf /Vj (middle panel), and (J/M) f /(J/M)i (bottom panel) is plotted against the final to initial 
disc mass ratio Mf /Mi for various combinations of nij and Zi. 



3.3 Dependence on baryonic matter distribution 

It is instructive here to additionally examine the dynamical re- 
sponse of an exponential disc for which the surface mass density 
distribution is given by 



S(r) = ^ exp(-r/r s ), (26) 

where r s is the scale length of the disc. This choice of E(r), which 
allows no analytic expression of W, is known to be more appropri- 
ate to spiral galaxies than the Kuzmin disc. 

For the purpose of comparison, we equate the effective ra- 
dius of the exponential disc to that of the Kuzmin disc, or r e = 
1.68r s = V3rt- Then, the dynamical response of the exponential 
disc embedded in the NFW halo is calculated by numerically inte- 
grating W, and the result is shown in Figure[5]in the same way as 
in Figure Q] 

It is evident from these figures that, given rrn and Zi, the dy- 
namical response of the exponential disc is only slightly stronger 
than that of the Kuzmin disc. Namely, the difference is very small 
for the impulsive mass loss and is much smaller for the adiabatic 
mass loss. This is understood, because the Kuzmin disc is more or 
less similar to the exponential disc except for the central region. 
Therefore, without loss of generality, we use the Kuzmin disc in 
Section 4 for the study of dynamical response in spiral galaxies 
over a full range of magnitude observed. 



4 SCALING RELATIONS MODIFIED BY DYNAMICAL 
RESPONSE 

In this section, using the Kuzmin disc, we investigate the effects of 
dynamical response in predicting the rotation velocity versus mag- 
nitude relation (TFR) and the disc size versus magnitude relation. 

We begin with the assumption that dark matter haloes 
in the CDM universe ob ey the scaling relations (e.g. 
iNavarro. Frenk & Whitdl 19971) : 

M d oc \Q rc oc r' d . (27) 

The baryonic matter in an individual halo falls in to form the galaxy 
disc and a part of such fallen gas goes back to the halo region ow- 
ing to the gas removal induced by SNe in the disc. Cosmologi- 
cal iV-body/SPH simulations have shown that along this scheme of 
galaxy formation the above scaling relations of dark matter haloes 
almost hold in brigh t spiral galaxies where the effect of SN feed- 
back is minimal (e.g. lSteinmetz & Navarrol 19991) . Consequently, in 
this paper, we assume that initial discs before the gas removal obey 
the same scaling relations as dark matter haloes, instead of enter- 
ing into detailed discussions of physical processes involved. This 
assumption accords with our use of m,: = Mi/Md and Zi = r^/r^ 
as constant parameters, yielding 

Mi oc V? oc rf, (28) 

where the subscript "i" stands for the quantities of initial discs be- 
fore the gas removal as in previous section. 

Next, we express the quantities of final discs after the gas re- 
moval, to which the subscript "/" is attached. Considering that disc 
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Figure 2. Dynamical response of the Kuzmin disc embedded in the NFW halo. Same as in Figure [T] but for the dependence of dynamical response on the 
concentration parameter c with m; = 0.2 and Zi = Tj/ (r d /c) = 0.05. 




Figure 3. Dynamical response of the Kuzmin disc embedded in the dark halo. Same as Figure[T] but for the truncated spherical halo with the homogeneous 
density distribution. 
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Figure 4. Dynamical response of the Kuzmin disc embedded in the dark halo. Same as Figure [T] but for the truncated spherical halo with the 1/r density 
distribution. 




log (M/Mj) log (Mj/Mj) 



Figure 5. Dynamical response of the exponential disc embedded in the NFW halo. The results are shown in the same way as in Figure[T] for direct comparison 
with the results for the Kuzmin disc. We note that the scale length r B of the exponential disc is related to rj in z via r s = (\Z3/l-68)r(,. 
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mass is a sum of remaining gas and stars M< = M gas + M, be- 
fore the gas removal and is left with stars Mf = M* after the gas 
removal, and introducing the strength parameter for SN feedback 
f3 = Mgas/M* for which [3 <C 1 for massive and normal galaxies 
and /3 ^> 1 for dwarf galaxies, we obtain 



M f 
M~ 



1 



1+/?' 



(29) 



Follow ing the definition in equation (10) of Nagashima et al. 
d2005h . we set 



What) 



(30) 



where Vhot is an effective SN-feedback strength in units of velocity, 
and ahot is a constant power index. 

Our analysis below confines to the NFW halo, because the ho- 
mogeneous and 1/r density distributions are of academic interest 
only. Values of the parameters used are Vhot = 150 km s -1 and 
cthot = 2 and 4 for the SN-feedback strength, while rm = 0.2 and 
Zi = 0.05, 0.1, and 0.2 for the initial disc. 



4.1 Tully-Fisher relation 

Theoretical TFR or the Vf — Mf relation is obtained by a set of the 
following equations: 

(31) 



Mf oc 
and 
Vf- 



1 + /3' 



Vi 

Vil[mi 



«; 1/(1 + 0)] 



(no response) 
(with response) , 



(32) 



where I (a, b; x) is given in equation l |16l >. 

This relation is compared with the data of rotation velocity 
Vrot and absolute mag nitude Mi of spiral galaxies in the I band 
taken from the table of Mathewson, Ford & Buchhornl dl992h . and 
the results are shown in Figure|6]for a^ot = 2 and in Figure|7]for 
cthot = 4, together with the data. In each of these figures the left 
and right panels are for the impulsive and adiabatic mass losses, 
respectively. 

Here, assuming a constant baryonic mass to I-band light ratio 
Mi jh\ and equating Vi to the observed rotation velocity V ro t for 
bright galaxies, we have set the power index 7 = 3 in equation 
d31t by adjusting the Vf — Mf relation along the horizontal axis in 
Figures[6]and[7]to fit to the observed TFR in the bright end. Our set- 
ting of 7 = 3 is consistent with the scali ng relation of dark haloes 
obtain ed by iV-body CDM simulations dNavarro. Frenk & Whitel 
as well as the scaling relat ion of spiral galaxies obtain ed by 
TV-body/SPH CDM simulations dSteinmetz & Navarrdl 19991) . 

We see from these figures that the V/ — Mf relations with 
no dynamical response and with no dark halo deviate significantly 
from the faint data, while the dynamical response with dark halo 
improves the fit to the data. In particular, the Vf — Mf relations 
that well agree with the data over the full range of Mi observed are 
those of Zi — 0.2 and cthot — 2 — 4 for the impulsive mass loss, 
and those of Zi = 0.1 — 0.2 and cthot = 2 as well as Zi = 0.2 and 
cthot = 4 for the adiabatic mass loss. 



4.2 Disk size versus magnitude relation 

Theoretical disc size versus magnitude relation or the 77 
relation is obtained by a set of the following equations: 



Mf oc 



and 



rf = 



1 + /3' 



Vi (no response) 

riH[irii, Zi] 1/(1 + /3)] (with response), 



(33) 



(34) 



M, 



where H (a, b; x) is given in equation dl5l >. Here, Vi(ri) is obtained 
from the relation of (Vi/Vhot) 1 = (ri/rhot) e according to equa- 
tion l |28l l, where Thot corresponding to Vhot is obtained from the 
rotation velocity versus magnitude relation and the disc size versus 
magnitude relation. So far we have used the characteristic radius n 
to represent the disc size. However, we hereafter use the effective 
radius r e = v3rt of the Kuzmin disc to be consistent with the 
usual observational definition. 

This relation is then compared with the data of effective disc 
radius r e and absolute mag nitude Mb of sp iral galaxies in the B 
band taken from the table o f llmpevetalJd 19961) . and the results are 
shown in Figure [8] for cthot = 2 and in Figure [9] for cthot = 4, 
together with the data. In each of these figures the left and right 
panels are for the impulsive and adiabatic mass losses, respectively. 

In order to avoid any systematics arising from our use of the 
velocity data in the I band and the size data in the B band, we first 
transform the theoretical relation from the I band to the B band 
by applying the empirical formula Mb = O.86A/1 — 1.31, then 
we compare it with the data of r e and Mb in the B band. Such 
a transformation formula is obtained b y eliminating Vrot fro m the 
observed TFRs in the B and I bands dPierce & Tullv|[l988h , and 
gives B— 1~ 1.2 for dwarf galaxies with Mi ~ —18 + 5 log h, and 
B— 1~ 1.6 for bright galaxies with Mi ~ —21 + 51og/i. These 
colors agree well with the observed colors (B— V) = 0.4 — 0.5 
and (V-I) = 0.7 - 0. 8 for dwarf ga laxies dSchombert et al J 19951 
1 19971 : lMakarova|[T999l : lO'Neil et alJ koOO). and with the synthetic 
colors of (B-V) = 0.5 - 0.6 and (V-I) = 1.0 - 1.1 f or bright 
Sb and Sc galaxies from Table 3a of Fukug ita et al.l dl995h . 

Here, applying the above transformation formula and allow- 
ing a small, constant offset of logr e (B) — logr e (I) between the 
effective radii in the B and I bands, we have set the power index 
e = 2 in equation i33\ by adjusting the 77 — M/ relation along the 
horizontal and vertical axes in Figures [8]and|9]to fit to the observed 
disc size versus magnitude relation in the bright end. Our setting 
of e = 2 is consistent with the scaling relation of dark haloes pre- 
dicted by a power-law spectrum of initial density fluctuations with 
an index n = —2 which is a reasonable approximation to t he CDM 
spectrum on galaxy scales (Navarro, Frenk & White 1997). 

We find that the 77 — Mf relation is much more sensitive to 
the parameters than the Vf — Mf relation, and many of the pre- 
dicted relations cannot reproduce the faint data at all. In particular, 
all the relations considered for the impulsive mass loss should be 
rejected by the data. On the other hand, for the adiabatic mass loss, 
the relations of Zi = 0.1 — 0.2 and cthot = 2 agree with the data 
over the full range of Mb observed, and the relation of z; = 0.2 
and cthot = 4 agrees with the data brighter than Mb = — 16 only. 



5 RESULTS 

Current SA models, where the mass loss from individual galax- 
ies by SN feedback is taken into account without considering dy- 
namical response, are known to overpredict the rotation velocity of 
dwarf spiral galaxies (thick solid lines in Figures[6]and[7} and at the 
same time underpredict their disc size (thick solid lines in Figures 
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Figure 6. Theoretical I-band Tully-Fisher relations with and without dynamical response of the Kuzmin disc embedded in the NFW halo. The results with 
a hot = 2 are shown by various thick lines in the left panel for the impulsive mass loss and in the right panel for the adiabatic m ass loss, while the result with no 
halo is shown by thin line only for the academic purpose of comparison. Filled circles represent the data taken from the table of Ma thewson. Ford & Buchhornl 
is the rotation velocity, M i is the absolute I-magnitude, and h is the Hubble constant Hq defined as h = Hq / 100 km s 1 Mpc 1 . 




Figure 7. Theoretical I-band Tully-Fisher relations with and without dynamical response of the Kuzmin disc embedded in the NFW halo. Same as in Figure 
[6] but for a hot = 4. 



[8] and [9}. We could easily imagine that such discrepancies are re- 
solved by the dynamical response of a virialized system associated 
with weakening of the gravitational potential, most likely owing to 
the mass loss that we consider in this paper. 

In fact, there exist several solutions with dynamical response 
that account for the observed relation between rotation velocity and 
absolute I-magnitude (Subsection 4.1) and the observed relation be- 
tween disc size and absolute B-magnitude (Subsection 4.2). How- 
ever, these two relations can only be reproduced simultaneously 
in a limited range of parameter space for the adiabatic mass loss, 



namely m.j = 0.2 and Zi = 0.1 — 0.2. For the SN feedback, the 
weak case of Qhot = 2 is preferable to the strong case of Qfhot = 4. 

Such an almost unique solution with Qhot = 2 for the adia- 
batic mass loss gives the final ratio of z/ — 0.3 — 0.4 for dwarf 
galaxies. This value of Zf agrees with observations in the local uni- 
verse (see Section.^. The value of (J/M) f/(J/M)i beyond unity 
suggests that the mass loss does not accompany the loss of angular 
momentum. 

The results above are obtained for dark haloes with a concen- 
tration parameter of c = 10, but it is known that c is predicted 
to vary from 10 for large haloes to significantly larger values of 
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Figure 8. Theoretical B-band disc size versus absolute magnitude relations with and without dynamical response of the Kuzmin disc embedded in the NFW 
halo. The results with a/, oi = 2 are shown by various thick lines in the left panel for the impulsive mass loss and in the right panel for the adiabatic mass 
loss, while the resu lt with no halo is shown by thin line only for the academic purpose of comparison. Filled circles represent the data taken from the table of 
llmpevetalJh996l) .rv is the effective disk radius, Mb is the absolute B-magnitude, and h is the Hubble constant Hq defined as h = Hq/100 km s 1 Mpc 1 . 




Figure 9. Theoretical B-band disc size versus absolute magnitude relations with and without dynamical response of the Kuzmin disc embedded in the NFW 
halo. Same as in Figure [8] but for 



20 — 30 for small haloes teullock et alj2001bl : lMaccio et alj2007h . 
It is then worth examining this dependence with other parameters 
fixed. The results by varying c for the Tully-Fisher relation and the 
disk size versus magnitude relation are shown in FiguresfTOlandl 1 1 1 
respectively. As seen from these figures, contrary to the case of im- 
pulsive mass loss, the variation in c does not make much difference 
in the results for the adiabatic mass loss, which validates our con- 
clusion in favor of adiabatic mass loss, irrespective of concentration 
of the mass included in dark haloes. 



6 DISCUSSION 

6.1 Losses of mass and angular momentum 

The angular momentum distribution in dark haloes is well 
established theoretically based on the CDM model (e.g. 
ICatelan & Theuns!ll996l : iNagashima & Gouddll998h . It is natural 
to assume that the baryonic components acquire the same amount 
of specific angular momentum as that of their host haloes, because 
the cosmic tidal field provides the ang ular mome ntum to dark mat- 
ter and baryons in the same manner ( Whit e] 1984h . This assumption 
has been a basis of most theories of galaxy formation. 
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Figure 10. Theoretical I-band Tully-Fisher relations with and without dynamical response of the Kuzmin disc embedded in the NFW halo. Same as in Figure 
[6] but for the dependence on the concentration parameter c. The results are shown by various thick lines in the left panel for the impulsive mass loss and in the 
right panel for the adiabatic mass loss. 




Figure 11. Theoretical B-band disc size versus absolute magnitude relations with and without dynamical response of the Kuzmin disc embedded in the NFW 
halo. Same as in Figure [5] but for the dependence on the concentration parameter c. The results are shown by various thick lines in the left panel for the 
impulsive mass loss and in the right panel for the adiabatic mass loss. 



From observations of dwarf spiral galaxies, 
Ivan den Bosch et all d2001ah found that their mean (or total) 
specific angular momenta are very similar to those of their host 
haloes predicted by those theories. Moreover they found that the 
the disc to halo mass ratio is smaller than the baryon fraction in 
the Universe. This means that only a part of baryons forms discs 
whereas it keeps an amount of specific angular momentum similar 
to that of host haloes. They also examined the specific angular mo- 
mentum distribution in individual haloes, and surp risingly found 
that it i s inconsistent with theoretical prediction by iBullock et al.l 
l l2001al) based on high-resolution iV-body simulations. 



Secondly, the observed surface density distribution of dwarf 
spiral galaxies is described as an exponential distribution (e.g. 
iMathewson. Ford & Buchhornll 1992b . different from that expected 
from the CDM model which provides much more centrally con- 
centrated discs when compared with the exponential distribution. 
This cannot be solved even if SN feedback is simply considered 
( Ivan den Bosch1l2001bh . 

One possible solution to the above discrepancies should be 
obtained by considering the mass loss from the central region of 
discs, where the gas should have the specific angular momentum 
lower than that in the outer region. SN feedback driving the mass 
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loss without the loss of angular momentum leads discs to increasing 
their specific angular momentum, i.e., (J/M) //( J/M)i > 1. This 
situation is consistent with that considered in this paper, as shown 
in the previous section. 

So far we have assumed rotation-dominated discs throughout 
gas removal. With this assumption, the specific angular momentum 
always increases after the mass loss, yielding (J/M) f / (J/M), > 
1. On the other hand, the transformation of rotation-dominated 
discs to velocity dispersion-dominated discs, which is the puff- 
up transformation, leads discs to decreasing the specific angular 
mom entum, yielding (J/M) f /(J/M), < 1 dBiermann & Shapiro! 
1979). This situation should occur when the gas in the outer region 
is removed where the specific angular momentum is large. Such 
a way of gas removal w ould be caused by the tidal stripping in 
clusters of galaxies (e.g. iGunn & iFuiita & Nagashima 

1 19991 ; lOkamoto & Nag ashima 20031). iBiermann & Shapiro! dl979t) 
proposed that this mechanism should be of lenticular galaxy for- 
mation. If so, the fraction of SO galaxies in clusters will be higher 
than in fields, and the rotation velocity of dwarf galaxies in clus- 
ters will be systematically lower than that in fields. This should be 
observationally clarified in future. 



6.2 Density profiles of dark haloes 

We have assumed three types of density profiles of dark haloes, 
such as the NFW profile with an inner slope of -1 and an outer 
slope of -3 (equation 1221 . and two power-law profiles with a slope 
of (homogeneous) and -1 (r~ ). 

The NFW profile is suggested by iV-body simulations, 
and many authors claim that the inner slope will become even 
steeper after the conde nsation of baryons owing to the cooling 
(Blumenthal et al. 198 ^), which is c a lled a s the adiabatic contrac- 
tion (but see ISellwood & McGau gh (2005) for the effects of ran- 
dom motions). In fact, this implicitly assumes that the cooling 
timescale is longer than the relaxation timescale for dark haloes 
to settle into the NFW profile, that is, baryons cool and shrink af- 
ter virialisation. On galaxy scales, however, this would not be the 
case. The cooling timescale is much shorter than the dynamical 
timescale, and the relaxation timescale is similar to or longer than 
the dynamical timescale. Thus it is not unnatural to assume that a 
mixture of dark matter and cooled baryons virializes. In this case, 
there would be galaxy discs within dark haloes with the NFW pro- 
file, without undergoing the adiabatic contraction. This might be an 
opposite limiting case to the adiabatic contraction, but recent X-ray 
observation of the intracluster medium h as found that there is n o 
evidence of contraction of the dark halo dZappacosta et alj |2006). 
Since the central cD galaxy is massive enough to suppress the SN 
feedback, we do not need to expect the expansion due to dynami- 
cal response to SN-induced gas removal. Thus it is possible to say 
that a simple assumption of adiabatic contraction does not work 
in reality. This should justify our use of the NFW profile for dark 
haloes that surround galaxy discs. To clarify this, of course, high- 
resolution hydrodynamical simulations includ ing the gas cooling 
processes are required (e.g. iGnedin et alj|2004l) . 



6.3 Hierarchical formation of galaxies 

In Section 4, our results of dynamical response to SN-induced 
gas removal within dark haloes are applied to scaling relations for 
galaxies, under the assumption that simple scalings among mass, 
velocity, and size of baryonic components have been set up before 



the gas removal. It is reasonable to consider that these are inferred 
from observations of scaling relations for massive galaxies where 
SN feedback is not effective owing to their deep gravitational po- 
tential wells. 

Thi s seems to be somewhat in contrast to the approach of, for 
example. lvan den Boscbl 1200 1b). in which the direction of angular 
momentum vector is assumed to be invariant. However, it has been 
shown that the direction can be moved by contiguous accretio n 
of dark matter dNagashima & Goudafl998llSaitoh & Wadal2004l) . 
Therefore, it is not assured that the specific angular momentum of 
discs is similar to that of host haloes averaged over whole regions. 

We thus believe that it is reasonable to use the scaling relations 
for massive galaxies as the initial state for the dynamical response 
to gas removal. 



7 SUMMARY 

We have analyzed the dynamical response to SN-induced gas re- 
moval on rotation velocity and disc size of spiral galaxies, explic- 
itly taking into account the underlying gravit ational potential wells 
made by dark matter. This is an extension of Nagashi ma & Yoshill 
(2003) in which spherical galaxies were considered, but for thin 
discs in spherical dark haloes to describe spiral galaxies realisti- 
cally. Similarly, expansion of discs decreases their rotation velocity 
and increases their size owing to the dynamical response to gas re- 
moval even within dark haloes. As shown in Section 3, the dynami- 
cal response provides unavoidable effects on the evolution of spiral 
galaxies, particularly on dwarf spiral galaxies because of their shal- 
low gravitational potential wells. 

We have examined the effects of dynamical response on scal- 
ing relations among mass, rotation velocity, and disc size of spi- 
ral galaxies. Since it is complicated to make such relations within 
a framework of hierarchical formation scenario of galaxies as 
discussed in Section 6, we use scaling relations which massive 
galaxies satisfy as the initial state before the gas removal. The 
SN-induced gas removal decreases the mass particularly of dwarf 
galaxies, thus distorting their scaling relations significantly. Such 
a distortion in the TFR as well as the disc size versus magnitude 



els of galaxy formation (e.g. 


Cole et al.l200d:lNagashima & Yoshiil 


12004 Nagashima et alj2005 


), in which the dynamical response on 



spiral galaxies is not taken into account. We have shown that the 
dynamical response considered in this paper is able to reproduce 
the scaling relations as observed. 

The hierarchical clustering scenario provides complicated 
galaxy formation processes. There are many physical processes 
such as radiative gas cooling, star formation and SN feedback, 
within merging dark haloes. Furthermore, while the gas removal 
occurs mostly due to SN feedback in massive galaxies, other pro- 
cesses including photoheating by the U V background may un- 
likely scale the same way in dwarf galaxies ( Nagashim a et al J 19991 : 
ISommervillel2002l ; lBenson et al. 2003; Ho eft et alj2006l) . In addi- 
tion, the formation history is also complicated. For example, discs 
accrete the gas which is once expelled from galaxies. It means that 
discs get bigger even after the dynamical response has exerted upon 
them. Then, luminous discs tend to cause large extinction by dust, 
which would bend the TFR. To investigate the observed properties 
of galaxies, therefore, we need to fully incorporate the effects of dy- 
namical response into realistic models of galaxy formation like our 
SA models. We will study the physical origin of scaling relations 
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including the TFR in more realistic situations in a forthcoming pa- 
per. 
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APPENDIX A: GRAVITATIONAL POTENTIAL ENERGY OF TWO-COMPONENT SYSTEM 

In this section we provide analytic expressions of f(z) and q(z) for the gravitational interaction potential energy of the Kuzmin disk of 
baryons embedded in various density distributions of dark matter halo. 



Al The NFW dark halo 

The NFW density distribution and the corresponding gravitational potential are given respectively by 

p(r) = p d c 3 
and 

<£>(r) = ~4tvG?^ In ( 1 + -V 
r V r d J 

Then, we obtain the total mass of dark halo 



cr ( _ cr 
— 1 + — 

r d V r d 



y 

r d ) 



M d = 4ivp d r d 



In (1 + c) 



1 + c 



and the analytic expression of the following function: 



f(z)=c 



In (1 + c) 



l + c 



1 , cz 1 
— In — H ; 

cz 2 czVc 2 z 2 + 1 



(VC 2 2 2 + 1 + 1){VC 2 Z 2 + 1 + CZ) 



(Al) 



(A2) 



(A3) 



(A4) 



For this case, we cannot obtain the exact analytic expression of q(z). Instead, we obtain an approximate expression of q(z) around z = as 
follows: 



q(z) ~ c 



In (1 + c) 



1+c 



/ 7 2 , cz\ 



4 2 3 



(A5) 



A2 The homogeneous dark halo 

The homogeneous density distribution and the corresponding gravitational potential are given respectively by 

p(r) = p d 0(r d - r), 
and 



$(r) = -4-rvGp d 



2 2 / 2 2 3 



Then, we obtain the total mass of dark halo 

4 3 
M d = -nr d p d , 

and the analytic expressions of the following functions: 

/m=^ 2 +!+-- V+i) 3/a , 

Z z z 

and 



q{z) = -y/z 2 + l(2z 2 - 1) + 2(z 3 - 2) + - sinh _1 (z). 

z 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



A3 The r 1 profile dark halo 

The 1/r density distribution and the corresponding gravitational potential are given respectively by 

p(r) = p d r -^Q(r d - r). 

and 



<E>(r) = -A-n:Gp d r d 



2r d 



+ 1 



1 + 



) 0(r - r d ) 



(All) 



(A 12) 



Then, we obtain the total mass of dark halo 



M d = 2-Kr dPd , 



(A13) 
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and the analytic expressions of the following functions: 

f{z) = 2 + z\nz- zln(l + ^/ z 2 + 1) + - - ^ * + 1 , (A14) 

z z 

and 

q(z) = -4 + | V^ + l + I" Sinn" 1 (*) + In - |« 2 ln(l + + (A15) 
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